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In aqueous solutions, dissolved ions interact strongly with the surrounding water, thereby mod- 
ifying the solution properties in an ion-specific manner. These ion-hydration interactions can be 
accounted for theoretically on a mean-field level by including phenomenological terms in the free 
energy that correspond to the most dominant ion-specific interactions. Minimizing this free energy 
leads to modified Poisson-Boltzmann equations with appropriate boundary conditions. Here, we 
review how this strategy has been used to predict some of the ways ion-specific effects can modify 
the forces acting within and between charged interfaces immersed in salt solutions. 



I. INTRODUCTION 

Resolving the interactions of ions in aqueous solu- 
tion can be a daunting task due to their complex inter- 
dependency. Consequently, the theory of electrolytes has 
been an active field of research for the last century and 
continues to be so even today. In the early 20" 1 cen- 
tury pioneering achievements by Gouy, Chapman, Debye, 
Hiickel and Langmuir resulted in the so-called Poisson- 
Boltzmann (PB) theory. The starting point of this mean- 
field theory was to consider point-like ions immersed in 
a solvent that was modeled as a continuum dielectric. 
Neglecting charge correlations and fluctuations, the elec- 
trostatic energy and the ionic entropy led to impressive 
predictions on the osmotic pressure, ion activities, and 
ion profiles close to charged interfaces and electrodes. 
These findings still form the basis of much of our cur- 
rent understanding of electrolyte solutions [lj. 

Even earlier, however, in the late 1800s Hofmeis- 
ter (among others) encountered ion-specific phenom- 
ena when he followed the effects of various salts on 
macromolecular interactions in solutions Q . Specifically, 
Hofmeistcr showed that certain monovalent ions (such as 
fluoride and chloride) are more effective at precipitat- 
ing proteins ("salting out") than others, such as bromide 
and iodide. This ion-specific capacity is often described 
in terms of the "Hofmeister series" and is correlated with 
other properties of aqueous ionic solutions, such as their 
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surface tension and solution viscosity. To explain these 
observations, one must consider additional ion-ion and 
ion-solvent interaction, which are not directly and exclu- 
sively related to the ion charge. To quote Hofmeister 
himself 0: "It also can be expected that the precipi- 
tating capacity of salts is parallel to other physical and 
chemical properties, if [...] these properties are depen- 
dent on the water adsorbing capacity of the salts" . 

In the last decades, the number of observations 
that found similar ion-specific interactions has steadily 
grown 0, H|. These effects are found to influence the 
interactions of surfactant micelles, lipid-bilayer mem- 
branes, proteins, DNA molecules, and more. As an exam- 
ple we mention charged and net-neutral lipid membranes 
that can swell or remain unswollen depending on the type 
of counterions @ . This observation can be related to the 
"stickiness" of larger and more polarizable ions at the 
vicinity of membranes. Ion-specific interactions result in 
a net exclusion or inclusion of ions at interfaces, which 
then can drive macromolecular interfaces, such as lipid 
membranes, to bind or separate. An extensive overview 
is beyond the scope of this review, but the interested 
reader should consult the recent surveys 0, Q ■ 

Beyond their polarizing effect on water molecules due 
to their charge, ions are themselves polarizable, and 
therefore interact with water and interfaces through van 
der Waals (dispersion) forces. These forces are ion spe- 
cific, and depend on the ion polarizability. In addition, 
ions are far from being simple point-like particles, and 
also show steric interactions that depend on ion size, as 
quantified by their effective ionic radii. Once we con- 
sider more complex ions that are composed of several 
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nuclei, interactions between permanent dipoles should be 
examined. Furthermore, these ions can have more com- 
plex structure and manifest additional forces, such as hy- 
drophobic interactions between water and nonpolar parts 
of the ions. Since all these interactions are mediated by 
the solvating water molecules, we will jointly refer to all 
these additional interactions under the broad name of 
hydration forces, without fully defining their microscopic 
origin. The focus of this review is on the ways hydration 
interactions can be incorporated and used through sys- 
tematic phenomenological modifications of the PB the- 
ory. 

The complexity of hydration interactions makes it al- 
most impossible to account for all their sources ab initio. 
The structure of liquid water and hydrogen-bonded net- 
work is in itself highly complex to model and predict. 
Even strategies that rely on molecular dynamics simula- 
tions inherently include some level of approximation [9(. 
As an alternative, we discuss a phenomenological strat- 
egy, where the free energy is our primary tool. The prin- 
cipal idea is that any interaction, including hydration 
interactions, can be added to the PB free energy, and 
its effect on the osmotic pressure and ion profiles can 
be subsequently assessed. Schematic drawing of some 
of these additional interactions is presented in Figure [1] 
In these endeavors we will, however, confine ourselves 
to short-range hydration interactions that allow for the 
most straightforward analysis. Similar strategies have 
been successfully used in other instances where details of 
specific interactions are not fully known, such as in lipid 
membranes and liquid crystals. 

We start in Sec.|TI]by presenting the PB free energy of 
ionic solutions, and the corresponding expressions of the 
density profile and osmotic pressure. Then, in Sees. IIIII 
and IIV1 we show how additional terms accounting for 
other non-electrostatic interactions lead to important ef- 
fects and have consequences to macromolecular interac- 
tions. 



II. GENERAL APPROACH 



We base our presentation on the PB free energy as the 
primary starting point fioj , with its well-known mean- 
field limitations For simplicity, let us focus on a 
simple model system composed of two charged planar 
surfaces (of infinite extent) located at z = ±D/2, im- 
mersed in an electrolyte bath with a dielectric constant e 
that contains several types of ions. Each surface carries 
a charge density of a (elementary charge per unit area) . 
The PB free energy J-pb of the ionic solution is composed 



of three parts: 



•Fpb - / d 3 r 
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d r at/j . 
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where e is the dielectric constant, So is the vacuum per- 
mittivity (in SI units), and ni and qi are, respectively, 
the number density and the valency of the i th ionic type. 
The first two terms correspond to the electrostatic en- 
ergy, while the sum of the third and fourth terms is equal 
to — Ts, where T is the temperature and s is the ideal 
entropy of mixing (per unit volume) of all mobile ions. 
Note that UiO? is the volume fraction of the i th species, 
where a is a microscopic length. The chemical potential 
Hi of the i th species, appearing in the fifth term, can be 
regarded as a Lagrange multiplier ensuring that the ions 
between the two surfaces are coupled with an ionic reser- 
voir of a given ionic strength. The last term reflects the 
surface contribution to the electrostatic energy. 

In order to facilitate the generalizations that will follow 
hereafter, we rewrite Eq. Q]in a more formal fashion: 



T= I & A rf v + J d'rf. 

d 3 r(w -Ts-J2 + / d Vs, ( 2 ) 

where fy = w — Ts — faiii is the generalized bulk 
free-energy density written as a sum of generalized en- 
thalpy and entropy densities, w = w({ni},ip' and 
s = s({ni}), respectively. In addition, f s contains all 
the surface contributions and depends on the surface de- 
grees of freedom such as the surface charge a and surface 
potential ip B . The decomposition of the free energy into 
a volume and surface parts assumes that the surface in- 
teractions are of short range. 

The Euler-Lagrange (EL) equations of the above free 
energy can then be decomposed into a volume and a sur- 
face contribution. The bulk equations are 

^_(^fv\_^fv =Q and ^V =0 (3) 
dz \ dip 1 J dip drii 

These EL equations reduce to the common PB equation 
in the case of IF = J-pp, as in Eq. Q] In addition, variation 
of the surface terms in Eqs. [T][2] results in 



dfs 
dtp 



se ip' s - 



(4) 
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Modifications in any of the free-energy terms in Eq. [2] 
lead to changes in the ionic density profiles and boundary 
conditions, and represent deviations from the standard 
PB form. 
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FIG. 1. Schematic representation of four possible additions to the Poisson-Boltzmann free energy that account for different ion-specific 
hydration interactions. Additional interactions can involve interactions between the ions themselves (top left), interactions of ions with surfaces 
and interfaces that are non-electrostatic (top right), as well as interactions with the solvent that itself could be either single (bottom right) or a 
multi-component mixture (bottom left). 



From the PB free energy, Eq. [TJ it follows that the 
pressure between the confining surfaces is 

i 

The pressure is composed of an attractive electrostatic 
term stemming from the Maxwell stress tensor (l2| . and 
a second term that resembles the van 't Hoff form and 
originates from the ideal entropy of mixing. 

The osmotic pressure is given by the difference between 
the pressure P at finite separation and its bulk value Pb 
(infinite separation). In the following we will refer to 
P as the osmotic pressure, omitting the reference to its 
bulk value. For systems with symmetric boundary con- 
ditions (equally-charged surfaces), the osmotic pressure 
is always repulsive [13} ■ Note that, in general, boundary 
terms in the free energy do not contribute explicitly to 
the pressure. 

Generalizations of the above electrostatic free energy 
are represented in Figure 1 and entail changes in all three 
of its parts, w, s and f s . Modifications in the enthalpy 
w are related to the dependence of dielectric response 
function e on solvent and solute compositions, as well as 



with other non-electrostatic interactions. The changes in 
the ionic entropy s are related with ionic steric interac- 
tions. Finally, modifications of the surface free-energy 
f a are wrought by the short-range non-electrostatic ion- 
surface interactions. The approach advocated here pro- 
vides a convenient way to classify non-electrostatic effects 
in terms of how they modify the standard mean-field free 
energy. 

The extended bulk free-energy fy can be written in a 
rather general but explicit form 

/y = -y£(K})(^) 2 +^eg i n i V + h({m}) , (6) 



where the dielectric response e ({n*}) takes a general form 
and depends on the ionic densities. The expression for 
fy includes the changes in the electrostatic interaction 
part (the first term), as well as in the third term, where 
h combines entropy and all other non-electrostatic inter- 
actions. 

The corresponding osmotic pressure can then be de- 
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rived as 



P = 



Sfv 
5D 



£o 
2 



Ede 



^ ^ ' Pin ■ 



drii 

2 I 

(7) 

Only modifications in fy enter explicitly into the pres- 
sure expression. However, surface modifications lead to 
changes in the density profiles, which in turn cause devi- 
ations in P. The function h ({rii}) encodes all the non- 
electrostatic couplings of the ionic densities. Hence, the 
osmotic pressure depends crucially on the assumed form 
of h as will be demonstrated in several cases hereafter. 



III. BULK: ION-ION AND SOLVATION 
INTERACTIONS 



In the following section we review several examples 
where the bulk free-energy, /y, is modified on the mean- 
field level, shedding new light on the notion of hydration 
interactions. 



A. Steric effects of hydration shells 

Why would one consider steric effects as a part of hy- 
dration interactions? In an aqueous solvent, the hardcore 
radius of ions is set principally not by ionic impenetra- 
bility but rather by the much larger and variable size of 
its hydration shell, composed of highly structured vici- 
nal water [5j. This structured shell excludes ions from 
the proximity of one another (top- left panel of Figure 1). 
The full account of this problem can be addressed by em- 
ploying liquid state theories, and requires extensive tools 
such as computer simulations or numerical solution of in- 
tegral equations [l4|. However, even with these elaborate 
tools, the state-of-the-art picture of solvent and solute 
structure is not well understood. 

An alternative approach is to use a more crude, yet 
simple, mean- field theory that accounts for steric effects 
within the Coulomb lattice-gas model. On this level the 
lattice constant is chosen to be equal to the hydration 
shell radius, a. The model accounts for density satura- 
tion close to charged objects and surfaces, but does not 
provide any insight into the microscopic packing struc- 
ture. 

The latti ce-g as model of Coulomb fluids has a venera- 
ble history |15| and has been apparently reinvented sev- 
eral times in the last century [161 ] , starting with Bikerman 
in 1942 [l7|. The lattice-gas entropy has the mean- field 
form 



TS({n^}) = - 



n,a 3 ln(a 3 rii) 



1 - ^ riia^j In ^1 - ^ n^a 3 ^ 




FIG. 2. Counterion density n as a function of distance z from 
a charged interface. Results of three models are presented: the 
standard PB theory (solid line); PB theory including ionic hy- 
dration shell, Eqs. [HE] of size a = 7.5 A (dashed line); PB 
theory including linearly varying dielectric constant, Eqs. 11011111 
with (3 = —20 M _1 (dotted line). The surface charge density is 
— e/50A~ 2 in all cases. 



and can be systematically derived via a field-theory ap- 
proach (18J. For simplicity, it is assumed that all ions 
have the same hydration shell size a but this condition 
can be relaxed [19| ■ The corresponding EL equations lead 
to a modified PB equation whose solutions show a pro- 
nounced saturation of the ion density in high potential 
regions, as can be seen in Figure[2] (dashed line) 20] . This 
limits the highest ionic density in the vicinity of charged 
surfaces, extending the electrostatic double-layer further 
into the bulk, as compared to the standard PB profile. 

An even older attempt to take into account steric ef- 
fects at charged surfaces was introduced by Stern in the 
1920s [21]. In Stern's model, ions are assumed to be 
fully excluded from a layer of a few Angstroms (the so- 
called Stern layer) adjacent to the surface, while outside 
this layer the profile is determined by the standard PB 
model. The lattice-gas and the Stern models are both re- 
lated to steric effects, although their ion-density profiles 
differ at the surface. 

The osmotic pressure for a Coulomb lattice-gas be- 
tween two charged surfaces was derived in Ref . |12| , hav- 
ing the form 



P = - 



f(^-^ln(l-E«.« 3 ) 



(9) 



(8) 



which can be reduced to the PB form of Eq.[S]in the limit 
of a small hydration shell (a 3 — > 0). 
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We note that steric exclusion, in general, augments 
the osmotic pressure between charged surfaces at small 
separation and provides a paradigm for the "hydration 
interaction" that stems from entropy contributions to the 
free energy. 



Ion-dependent dielectric response 

The formation of hydration shells around ions is driven 
by the strong local electrostatic fields [5|. Besides the 
steric effects (discussed in the previous section) due to 
the formation of the hydration shell, the dielectric re- 
sponse of the solution is also modified by the presence 
of ions |22j. Being dipolar, water molecules are strongly 
attracted to the ions and aligned themselves along the 
electrostatic field lines that originate from the ions. Con- 
sequently, the orientational degrees of freedom of water 
molecules are substantially reduced, leading to regions 
with lower dielectric response compared to pure water 
(bottom-right panel of Figure 1). This results in a decre- 
ment of the macroscopic dielectric constant as has been 
observed experimentally already in the late 1940's [23| . 
Furthermore, when the ions are distributed inhomoge- 
neously (for example, close to charged interfaces), the 
dielectric response will depend locally on the ionic pro- 
file. 

In the framework of a phenomenological mean-field 
free-energy, this effect [24[ can be accounted for by as- 
suming that the dielectric response depends linearly on 
the ionic density, n: 

e[n(z)] =s w +pn{z), (10) 

where e w is the dielectric constant of pure water, and f3 
is defined as a phenomenological coefficient. The linear 
relation for macroscopic dielectric constant was observed 
for several monovalent salts and for a fairly wide range of 
concentrations, see Ref. [24j and references therein. Note 
that only one ionic species (counterion only) is considered 
here for the sake of simplicity, but this can be easily 
generalized for any number of ionic species. 

The counterions now play a "dual" electrostatic role. 
On the one hand, they are attracted to the charged in- 
terface due to their charge, whereas, on the other hand, 
they are repelled due to their lower dielectric response 
relative to the pure solvent. The interplay between these 
two contributions leads to a decreased ionic density in 
the vicinity of the charged interface, when compared to 
standard PB theory. Specific results can be found in 
Ref. [13. 

The above mentioned dielectric response, Eq.[TUl leads 
to an explicit term, —(3£ n(ip') 2 /2, to be inserted in the 
bulk free-energy, fy of Eq. [6l In the prevalent case for 
monovalent salts, the linear coefficient is negative, /3 < 
0, and the free energy is increased due to the presence 
of the ions. The corresponding osmotic pressure can be 



obtained from Eq. [7] and is of the form 

P = ~ [e w + 2/3n] (V/) 2 + k B Tn . (11) 

The 2j3n (ip') 2 term together with the modified ionic pro- 
file (exclusion effects) results in an increase of the repul- 
sion between like-charged surfaces, as compared to the 
PB pressure, Eq. [5j as can be seen in Figure [2] (dotted 
line). This may appear as counterintuitive, because ex- 
clusion often entails effective attractive interactions be- 
tween interfaces (so-called depletion interactions). How- 
ever, it should be realized that the counterions should 
remain in the intervening space between equally charged 
surfaces, and they can only do so at a larger free-energy 
cost as compared to the standard PB theory. 

Ions in a dipolar solvent 

An alternative approach to the effect of ion-solvent 
interaction is to model water molecules as explicit, 
constant- magnitude point-like dipoles [25|, l26| . This adds 
another degree of freedom to the free energy of Eq. 
Since the dipolar degree of freedom can be integrated an- 
alytically, one obtains a free-energy ansatz along the lines 
of Eq.[5J with an additional polarization-energy term. In- 
stead of treating the solvent as a homogeneous medium 
with a macroscopic dielectric constant, the dielectric con- 
stant is calculated to be a function of the local electro- 
static field, ip': 

£cff = e + °-^Q{u) , (12) 

ip' 

where the bulk concentration of dipoles is Cd, Po is their 
dipolar moment, and u = poip' /k B T is a dimensionless 
energy variable accounting for the energy of a dipole po 
in an external electric field — ip'. The function Q{u) = 
coshu/u — sinhu/u 2 is related to the first order spherical 
Bessel function. The osmotic pressure for this model can 
be computed explicitly and leads to the form 

P = -^V) 2 + 

i 

sinhu 

-k B Tc d ug{u) + k B Tc d , (13) 

u 

where the third term is the dipolar electrostatic contri- 
bution, and the fourth term is the pressure of an ideal 
gas of dipolar particles with a concentration CdSmhu/u. 
The interplay between these two terms as well as the 
change in the ionic profile leads to a deviation from the 
PB pressure, making it smaller in the case of two oppo- 
sitely charged surfaces (25j . 

Ions in hydrogen-bonded aqueous solvent 

Another ion-specific effect explicitly considers the 
hydrogen-bond network of water molecules and its per- 
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turbation by the dissolved ions and confining sur- 
faces [27j. These types of considerations lead to finite- 
ranged hydration interactions and complicated couplings 
between the chemical structure of the confining sur- 
faces [28]. Of all the effects considered here these are 
the most difficult to describe theoretically in a consistent 
way, since new types of order parameters, in addition to 
ionic density and the corresponding electrostatic fields, 
need to be taken into account. Various types of phe- 
nomenological theories of finite-range hydration effects 
are mostly related to the microscopic Onsager-Dupuis 
model of ice [29| and are able to account for some facets 
of the structural correlations of the aqueous solvent as 
a consequence of its pronounced hydrogen-bonded struc- 
ture. 



Ions in binary mixtures 

Models accounting for ion-solvent interactions can also 
be generalized for the case of binary solvent mixtures. 
These are complex solvating environments, because the 
solution may respond to a solvated ion by varying the 
local solvent composition (bottom-left panel of Figure 
1). For example, a miscible mixture of high dielectric 
solvent (water) with a lower dielectric one (such as al- 
cohol) can respond to a solvated ion by drawing more 
water molecules around the ion on the expense of alcohol 
molecules in order to optimize the dielectric response [30| . 
Moreover, additional non-electrostatic interactions may 
preferentially drive one solvent to make closer contact 
with the ions. This general tendency is characterized by 
so-called preferential interactions, which compete with 
the entropic cost (or more generally, the free-energy cost) 
of concentrating one solvent around the ion. These inter- 
actions can, for example, be quantified by the differences 
in solvation energies of salts by different pure solvents and 
can even lead to macroscopic phase-separation of the liq- 
uids ("salting out") [3l|. Similar considerations can be 
applied to charged macromolecular interfaces. For ex- 
ample, recent experiments show that preferential water 
inclusion (preferential hydration) can lead to subtle, yet 
measurable changes in the forces acting between macro- 
molecules such as DNA [32j. 

At the simplest level, the local dielectric response in an 
A/B binary mixture depends linearly on the local solvent 
composition 4>a and 4>b as 



e(z) = 4>a{z)ea + 4>B{z)e B , 



(14) 



where due to solvent incompressibility (while neglecting 
the ionic volume fraction) , we require <f>A + <Pb = 1 ■ This 
form is supported also by experimental evidence [3l|. 

The total free energy is composed of three terms: 
T = -Fpb + ^"mix + -^soi- The first .Fpb term is the regu- 
lar PB expression, Eq. [1] with the dielectric constant e of 
the solvent substituted by the local dielectric response, 
e(z). The second term .F m i x is directly derived from reg- 
ular solution theory, and describes the ideal mixing and 



enthalpy of the solvent on the mean-field level in terms 
of the relative solvent composition <fi = <pB'- 



)ln0+ (1 - <j)) ln(l - <j>)] d 3 r 



k B T 



X 0(l-0)d 3 r, (15) 



where \ is the dimensionless interaction parameter 
(rescaled by ksT). The origin of the third term, Tso\, 
is the preferential (non-electrostatic) ion-solvent interac- 
tion, assumed to correspond to a bilinear coupling be- 
tween the ion densities, n^, and the relative solvent com- 
position </), so that 



ol 



aiTii4> d r 



(16) 



where oli is the preferential interaction parameter of 
species i. 

Minimizing the free energy with respect to <p, rn and 
ip y ields the corresponding equilibrium (EL) equations 
[33L l34l ] . Solving these equations results in ionic and sol- 
vent density profiles close to the charged surfaces. The 
osmotic pressure can be derived from Eq. [7] as a function 
of the inter-plate separation D and the experimentally 
determined by the bulk salt concentration rib, the bulk 
solvent composition and o^. It has the form 

+ fc BT (y^ m - ^3 ln (1 - ^ + 



k B T ^2 a i n i4> - ^<t 



(17) 



In the absence of preferential solvation, on = 0, and only 
as long as <C 1 and Ea > £b j small changes in osmotic 
pressure are expected, because density modulations of 
the solvent are small. However, considerable changes in 
the osmotic pressure are expected when preferential in- 
teractions are included [!§]; the effect is expected to be 
particularly large for separation between charged inter- 
acting surfaces of less than ~ 10 — 20 A. 



Solutions of antagonist ions 

Along similar lines as presented in the previous section, 
an interesting effect emerges by considering the peculiar 
case where co-ions and counterions bear an antagonist 
preferential solvation (3f| [36j|. This effect is of consid- 
erable importance to the properties of neutral interfaces 
where the two ionic species accumulate on different sides 
of an interface between two immiscible dielectric solvents. 
For example, we mention the oil/ water interface, where 
the hydrophobic ions accumulate on the oil side of the 
interface, and the hydrophilic ones on the aqueous side. 
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Under these conditions the surface tension is found to 
depend non-monotonously on the ionic bulk density, re- 
sulting from a negative contributing term proportional to 
^/rib, and a positive term linear in rn>. Note that in order 
to obtain this behavior even on the mean-field level, one 
must take into account image charge interactions, which 
become important near non-charged interfaces separat- 
ing two different dielectric solutions. Here these interac- 
tions were considered in a generalized Onsager-Samaras 
fashion and are further discussed in Sec. IIVI More re- 
cently 37], the implication of this effect on the phase 
diagram of aqueous mixtures was studied. It was found 
that due to this antagonist ion effect, hydrophilic ions 
tend to enhance phase separation between the two im- 
miscible solvents. 



IV. DIRECT IONS-SURFACE INTERACTIONS 

So far, we have discussed the contributions to the bulk 
energy of solvated ions, fy. Similarly, specific inter- 
actions can also occur at interfaces when the adsorbed 
ions interact favorably with the surface constituents (top- 
right panel of Figure 1). The main difference is in the 
lower dimensionality of the adsorbing surface. Below we 
review some examples of short-range surface-specific in- 
teractions that can be described via a modified surface 
free-energy, f s , (second term in Eq. [2]) , and examine the 
consequences of these modifications. 



Charge regulation 

In their pioneering model of charge regulation, Ninham 
and Parsegian 38] introduced the concept of a variable 
surface charge density that self-adjusts according to other 
system properties. The charged surface is not described 
by a constant charge density. Instead, the amount of 
charge is regulated according to association and disso- 
ciation of surface ionic groups. Considering such pro- 
cesses at thermal equilibrium yields a charge regulated 
boundary condition. The main addition in the Ninham- 
Parsegian model with respect to the regular PB model is 
the novel boundary condition, which is determined sclf- 
consistently and, in turn, fixes the amount of dissociated 
charge groups on the surface as a function of, e.g., the 
separation between the surfaces. This model proved to 
be very useful and was later applied in a plethora of other 
contexts [39| . 

Though originally the charge regulation was derived by 
requiring that the surface charging is at equilibrium with 
the bulk through a specified equilibrium constant for ion- 
surface binding 38] , it can be introduced equivalently by 
using a free-energy approach. Consider a surface com- 
posed of associated (neutral) and dissociated (charged) 
groups modeled as a two-component mixture. The free 
energy, f s (r), i(j s ) — at/j s + f cn t, is a sum of an electrostatic 



term and the surface entropy of mixing 



fe 



[77 11177 + (1 - 77) In (1 - 77)] , (18) 



where 77 is the area fraction of the dissociated (charged) 
groups at the surface, and a — erf /a 2 is the surface charge 
density. The boundary condition, Eq. |H is obtained from 
a variation of the free energy with respect to ip s , while 
minimization with respect to rj corresponds to the Lang- 
muir adsorption isotherm [4(| . Combining these two con- 
ditions leads to the Ninham-Parsegian charge regulation 
condition 1381: 



n 



l-r, 



exp 



(19) 



Since the boundary terms do not enter the pressure 
equation, Eq. the modification in the interactions be- 
tween two bounding surfaces is wrought only through the 
changes in the ionic profiles and the electrostatic field. 
General forms of surface free energies may lead to sur- 



face phase separations as discussed in Ref. 41] . As shown 



in the next section, such phase separation can couple to 
the spacing between the charged layers at equilibrium. 



Lamellar-lamellar phase transition 

One extension of the Ninham-Parsegian model consid- 
ers short-range non-electrostatic interactions between as- 
sociated and dissociated charged groups at the surface, 
and can be modeled by additional terms to f s [40] . This 
generalization is motivated by experiments by Zemb et 
al. [421 ] . where a lamellar-lamellar phase transition was 
observed in certain bilayer-forming lipids and surfactants, 
such as DDA (didodecyldimethylammonium halides), for 
a series of three homologous halides counterions: Cl _ , 
Br~ and I - . A discontinuous transition was found for the 
inter-lamellar spacing D as a function of applied osmotic 
stress, but only when counterions such as bromide were 
used. In contrast, for chloride there is no phase transi- 
tion and the isotherm, P(D), follows the PB prediction, 
while for iodide the lamellar phase did not disperse at all 
in water. 

There is a probable link between the surface activity of 
ions and the lamellar-lamellar phase transition character- 
ized by the discontinuous jump (for certain counterions) 
in the inter-bilayer separation. To suggest a possible ex- 
planation for these experiments, a free-energy approach 
can be used, with a surface term f s of the form 



fs = ai> s + f e 



1 



a s rj 



:Xsil 



(20) 



The modification in f s , as compared with the previous 
section, is included in the last two terms in the squared 
brackets. A non-electrostatic ion binding to the surface 
is modeled by the linear a s rj term, while the interaction 
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between the bound and dissociated groups is represented 
by the quadratic \xs r i L term. 

Minimizing the complete free energy leads to the regu- 
lar PB equation inside the bulk solution and to a bound- 
ary condition in the form of the Langmuir-Frumkin-Davis 
adsorption isotherm [4j|. The solution of the PB equa- 
tion with this boundary condition completely determines 
the counterion density profile rij(z) and electrostatic po- 
tential, and yields the osmotic pressure P of the form Eq. 
[5] as a function of separation D. 

The analysis of the solution shows that if the parameter 
Xs is large enough (typically s» lO/c^T), an in-plane phase 
transition can occur, with a coexistence region in the 
phase diagram between ion-adsorbed and ion-depleted 
phases. This transition, in turn, is coupled to the bulk 
phase transition causing a jump in the inter-lamellar 
spacing D as the osmotic pressure is changed. Quali- 
tatively, the behavior of this lamellar system, seen for 
different ions, can be understood in terms of the spe- 
cific counterion interaction with the charged surfactant 
bilayer. The Cl _ counterion always dissociates from the 
bilayer- forming DDA + surfactant, resulting in PB-likc 
behavior, and a continuous P{D) isotherm. For Br - , 
the dissociation is partial, leading to a first-order transi- 
tion and a coexistence between the two lamellar phases of 
two different D spacings. Finally, for the I - counterion, 
the ion always stays associated with the DDA+ surfac- 
tant and there is not sufficient repulsive interaction that 
can stabilize the swelling of the stack under any osmotic 
pressure. 

Bilayer curvature effects 

So far we have discussed rigid and flat surfaces. But 
often, interfaces are composed of lipids or surfactants 
that can be deformed at some free-energy cost. The 
adsorption process can then lead to a surface restruc- 
turing and deformation at equilibrium. For a multi- 
component layer, the compositional and elastic degrees 
of freedom are in general coupled, leading to a com- 
plex variational problem. Ion specificity is then man- 
ifested through charge regulation, through interactions 
between mobile associated and dissociated surface sites, 
and through other surface mediated interactions, such 
as bilayer elasticity. This is particulary relevant when 
considering the adsorption of large macroions onto oppo- 
sitely charged layers. 

Phenomenological terms added to the free energy are 
useful here too. The free energy can be written as a 
sum f a = / ent + / e i s . The additional term corresponds to 
an elastic energy in its Helfrich form, using the bending 
modulus k of the layer, 

/els = \ K i c ~ c o(v)} 2 (21) 

where c denotes the layer curvature, and c (rj) is the 
spontaneous curvature, depending on the local fraction 



of dissociated surfaces sites, r\. 

Using these modified equations, it has been possible to 
predict how charged macromolecules, such as DNA, pro- 
teins, and even viruses, can reshape oppositely charged 
membranes, and sometimes induce a macroscopic phase 
separation from a flat lamellar phase to a deformed 
one [i3 |. 



Surface tension of electrolyte solutions 

Another approach to study ion-surface interaction is 
to consider a non-electrostatic external field that has a 
direct dependence on the distance from the surface, and 
add this term "by hand" to the free energy, ending up 
with modified ionic profiles. The main difference between 
this method and those mentioned in the previous sections 
is the assumption of a certain functional form of the non- 
electrostatic external field, which is independent of other 
system properties. One famous example to this method 
is the Onsager-Samaras theory for the surface tension of 
electrolytes [45|. In this seminal work it was assumed 
that the image-charge potential can be regarded as an 
external potential that scales as ip- lm ~ exp(— 2X^ l z)/z 
where z is the distance from the surface and Ad is the 
Debye screening length. Taking into account V'im leads 
to a depletion layer in the vicinity of neutral interfaces. 

In recent works [46l - |48l |. a model was proposed to ex- 
plain the adsorption to an air/water interface of "soft" 
and large anions, such as iodide and bromide, and the 
depletion of "hard" and small cations, such as sodium. 
The model was based on two additional and independent 
external fields, absent in the standard PB theory. The 
contribution of the first stems from the cost of the elec- 
trostatic energy of a polarizable ion that is modeled as 
a sphere, positioned at the interface and immersed only 
partially in the liquid phase. This energy is found to be 
substantially reduced for a large soft anion. The second 
additional ingredient is the energy associated with ex- 
changing an ion with a water molecule in the hydrogen- 
bond network. Since the water molecule is energetically 
favored in the bulk, this will lead to an effective attrac- 
tion of the ion to the interface, which scales with the 
partial volume of the ion in the air phase. 

Taking these two independent external fields into ac- 
count leads to a dependence of the profiles on the ionic 
size and polarizability. The predictions suggest that the 
ion-specific surface tension agrees with the Hofmeister 
series 

In another work it was suggested that the specific 
ion property of making or breaking the water hydrogen- 
bond network can be modeled by an attractive poten- 
tial (square well) or a repulsive potential (square bar- 
rier) [49, H3| • Namely, the energetic favor or disfavor of 
ions to be part of the hydrogen-bond network leads to 
an effective attractive or repulsive interaction of the ions 
to the surface, respectively. Adding these external fields 
changes the ionic profile and induces an ion-specific de- 
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pletion or adsorption, depending on the choice of the ex- 
ternal potential. This non-monotonic behavior was also 
observed in surface tension experiments (5lj |. 

It is still questionable whether the simplified form of 
these external potentials can be fully justified. Further- 
more, in all models where independent external poten- 
tials are "inserted by hand" into the free energy, it is 
assumed that the additional interactions are additive. It 
remains an open question whether all these interactions 
are not intimately related to one another in a self consis- 
tent manner. 



Dispersion forces close to dielectric discontinuities 

Electrostatic interactions mediated by aqueous solvent 
are contingent upon the dielectric response function that 
encodes all the relevant structural properties of the sol- 
vent [52j |. The temporal dependence is always non-local 
and can be converted into an appropriate frequency- 
dependent dielectric response function, e(uj). The non- 
local spatial dependence, on the other hand, describes 
the orientational correlations between water molecules in 
a hydrogen-bonded network. For an infinite system it can 
be codified by a wave-vector dependent dielectric func- 
tion [53( | . However, additional approximations are needed 
in order to describe correlations in a system that is bound 
by interfaces and contains mobile counterions 

The link between the frequency-dependent e(ui) and 
the van der Waals (vdW) interactions is provided by the 
Lifshitz theory (55|. The dielectric models for the fre- 
quency response of aqueous solvents are well worked out 
and lead to quantitative estimates of these interactions 
for various geometries [5f|. Ninham and Mahanty were 
the first to suggest that the vdW interactions supply, at 
least in part, an explanation to the ionic specificity [57]. 
In particular, they suggested that the vdW interactions 
partially explain the origin of the water-structure contri- 
bution to the interactions between ions and surfaces. For 
ions interacting with surfaces the vdW potential in the 
non-retarded limit is given by ip v dw( z ) ~ B / z 3 where 
z is the distance from the surface, and the coefficient B 
depends on the polarizability of the ions. This sugges- 
tion has been picked up recently in several attempts to 
connect the Hofmeister series and other ion-specific ef- 
fects with the vdW interactions (see Refs. 0, [58| and 
references therein). 

At the most basic level the vdW interactions contribute 
an additional term to the mean-field potential leading 
to modified ionic profiles compared to the PB ones [59| . 
A more advanced mean-field continuum theory should 
take into account the dependence of the dielectric re- 
sponse on the local density of the ions [2J]. Thus, the 
frequency-dependent dielectric response could be writ- 
ten as s(z,lj) = e(ni(z) 7 uj). Note that the ion density 
dependence of the dielectric function has already been 
discussed in the previous sections, but the evaluation of 
the vdW interactions with a non-homogeneous dielectric 



profile is, in general, difficult (60j . 

An additional complication is due to the fact that not 
only the vdW interactions depend on rii(z), but the con- 
verse is true as well — the ionic profiles, rii(z), depend on 
vdW interactions. This coupled problem is even more dif- 
ficult to solve and, in general, may lead to non-monotonic 
density profiles [6l|. 

The vdW interactions add ionic specificity to the total 
free energy since e{m{z),u) depends crucially on the ion 
species. It is important, though, to note that this is only 
one part of the ion-specific effect. The other part has to 
be implemented within PB theory itself by consistently 
assuming that the static dielectric response depends on 
the ionic densities too, as advocated in previous sections 
of this review. Only taken together these two effects can 
constitute a consistent formulation. 



V. CONCLUDING REMARKS 

Ion-specific effects arc present in a broad range of 
chemical and biological systems. These effects are, how- 
ever, difficult to model since they usually involve com- 
plicated couplings that emerge from a large number of 
species and interactions. Here we have explored several 
models with a common theme — all are phenomenologi- 
cal and highlight one or two aspects of ion-specific inter- 
actions. Furthermore, all the models are discussed within 
a common theoretical framework on the mean-field level, 
leading to relatively simplified and intuitive picture of 
ion-specific effects. The drawback of this approach is the 
lack of accuracy and microscopic description that can be 
provided by other methods, such as simulations and de- 
tailed liquid-state theories. 

The various models can be grouped into two main 
types of non-electrostatic interactions: bulk and sur- 
face ones. For the first type we reviewed three non- 
electrostatic interactions. 

(i) The finite-size of ions, when non- negligible, leads 
to steric repulsion and to an upper bound on the ionic 
density that accumulates in the vicinity of charged sur- 
faces. As a consequence, the width of the double-layer 
is increased for larger ions and the osmotic pressure in 
between two surfaces, separated by a few nanometers, 
grows. 

(ii) Ions influence the dielectric response in their vicin- 
ity. Mostly, ions apply strong electrostatic field on neigh- 
boring water molecules, leading to a local decrement of 
the dielectric response. This, in turn, leads to an addi- 
tional effective interaction of ions with charged surfaces. 
The effect on the ion density profile is somewhat similar 
to that of large ions in the steric interaction model, and 
the double layer width is increased. Ion specificity in this 
model emerges via the strength of the ionic influence on 
the local dielectric response. 

(Hi) When the solvent is composed of several compo- 
nents, additional interactions lead to modifications of the 
ionic and solvent density profiles. In the model described 
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above, the change of local dielectric response due to lo- 
cal changes of solvent composition, and the difference in 
the solvation energy of the ions in different solvents lead 
to a decrease in the osmotic pressure as compared to a 
homogeneous solvent composition. 

The second type of non-electrostatic interactions are 
surface interactions and several cases were discussed. 

(i) Surface charge can be self-adjusted by various 
mechanisms, such as ionization and ionic dissociation. 
Such self-adjusted surfaces interact with the ionic solu- 
tion in a self-consistent manner, which leads to changes 
in the ionic density profiles and osmotic pressure. 

(ii) Adding an in-plane non-electrostatic interaction of 
the surface species to the charge regulation mechanism 
may further lead to a phase transition as a function of the 
osmotic pressure between two charged surfaces. Further- 
more, when coupling the charge regulation mechanism 
to the membrane curvature degrees of freedoms, struc- 
tural changes in the membrane can be observed when it 
is brought in contact with an ionic solution. 

(Hi) Considering long-range external surface interac- 
tions such as dispersion interaction, that result from the 



dielectric inhomogeneities in the system, leads to further 
modifications in ionic density profiles and osmotic pres- 
sure. Here too, the ionic specificity emerges from the 
ion-dependent dielectric response of the solution. 

Resolving a complete picture of ion-specific effects re- 
mains an ever-challenging task. In order to accomplish 
this task, it will be helpful to obtain first a better under- 
standing of each of the myriad experimental observations 
in terms of models based on simple ingredients. Then, 
combining the separate pieces of knowledge gained to- 
gether into a single and complete picture of ion-specific 
interactions may hopefully be achieved. 
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